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ABSTRACT 

The eclipsing binary system 2M 1938+4603 consists of a pulsating hot subdwarf B star and a cool 
M dwarf companion in an effectively circular three-hour orbit. The light curve shows both primary 
and secondary eclipses, along with a strong reflection effect from the cool companion. Here we present 
constraints on the component masses and eccentricity derived from the R0mer delay of the secondary 
eclipse. Using six months of publicly- available Kepler photometry obtained in Short Cadence mode, 
we fit model profiles to the primary and secondary eclipses to measure their centroid values. We find 
that the secondary eclipse arrives on average 2.06 ± 0.12 s after the midpoint between primary eclipses. 
Under the assumption of a circular orbit, we calculate from this time delay a mass ratio of q = 0.2691 
± 0.0018 and individual masses of M^d = 0.372 ± 0.024 Mq and Mc = 0.1002 ± 0.0065 Mq for the 
sdB and M dwarf, respectively. These results differ slightly from those of a previously-published light 
curve modeling solution; this difference, however, may be reconciled with a very small eccentricity, 
ecosw ~ 0.00004. We also report an orbital period decrease of P = (-1.23 ± 0.07) x 10~^". 
Subject headings: binaries: close — binaries: eclipsing — stars: individual (2M 1938+4603) — subd- 
warfs — techniques: photometric 



1. INTRODUCTION 

Around ten eclipsing hot s ubdwarf B (sdB ) star bina- 
ries are known presently fe.g.. lFor et al.|[2"010D : they have 
orbital periods from 2-4 hours, and the companions are 
typically low-mass M dwarfs. Although these subdwarfs 
outshine their companions at all visible (and most in- 
frared) wavelengths, the high albedo of the M dwarf and 
its large subtended solid angle as seen from the sdB star 
lead to a strong reflection effect, and, consequently, these 
systems exhibit both primary eclipses (when the M dwarf 
transits the sdB) and secondary eclipses (when the sdB 
blocks reflected light from the M d warf). Binary popula- 
tion synthesis (BI^S ) models (e.g.. iHan et"alll2002l [200l 
iClausen et al.|[20T^ show they are likely a product of a 
past common-envelope (CE) stage during which the sdB 
progenitor flUed its Roche lobe while on the red giant 
branch and lost most of its outer H envelope. Assessing 
the distribution of their orbital parameters and masses 
can test binary evolution scenarios and help constrain the 
parameterizations in BPS codes. Masses have been de- 
rived from light-curve modeling solutions for t he major- 
ity of known detached sdB+dM systems (e.g., iFor et aH 
mmh . but the accuracy of these results has not been 
tested thoroughly using independent techniques. 

Precise measurements of the eclipse timings in 
sdB+dM binaries provide an opportunity to measure 
the component masses in a relatively model-independent 
way. An observer watching one of these systems face-on 
from a great distance would see a syzygy of the sdB, its 
companion, and Earth every half orbital period, if the 
orbit is circular. An observer on Earth, however, has a 
different experience. Since the mass ratio does not equal 
unity, the primary and secondary stars will occult the 
light from the other star while at different displacements 
from the binary center of mass. That is, the eclipses 
as seen from Earth emanate from different line-of-sight 



distances. The secondary eclipses will arrive later than 
the halfway point between primary eclipses. (Here we 
use "primary eclipse" to refer to the eclipse of the pri- 
mary star, assumed to be brighter and more massive.) 
Using measurements of the orbital period, the sdB veloc- 
ity semi-amplitude, the inclination angle, and the time 
delay (R0mer delay) of the secondary eclipse, one can 
determine both masses, as if the system were a double- 
lined spectroscopic binary ([Kaplanl 120101) . The method 
is applicable to eccentric systems, too, but the orbital ge- 
ometry must be determined first since eccentricity shifts 
the relative timing of the two eclipses. 

2M 1938+4603, or KIC 9472174 (hereafter, 2M 193^1), 
is a bright (5= 11. 96) eclipsing sdB+dM binary that lies 
in the Kepler field and a suitable test case for the tech- 
nique described above. The system exhibits grazing pri- 
mary and secondary eclipses along with a strong reflec- 
tion effect from the cool secondary (jOstensen et al.ll2010L 
hereafter 010); the period is 3.0 h, and the orbit appears 
to be nearly circular. The analysis described by Kapla3 
(|2010f ) predicts the secondary eclipse should occur 2.35 
± 0.10 s after the midpoint between primary eclipses, as- 
suming the mass ratio from 01O's light curve modeling 
solution, g=0.244 ± 0.008. This time delay should be 
observable given sufficient measurements. After phase- 
folding a nine day long light curve obtained from Kepler 
during quarter (QO) of operations, 010 found evidence 
for more than 55 pulsational frequencies covering both 
the p- and g-mode domains in frequency space. Although 
the amplitudes are quite small (< 0.1 %), the rich pul- 
sation spectrum might permit asteroseismic models to 
establish the mass of the sdB. Thus 2M 1938 offers an 
opportunity to compare mass determinations for an sdB 
star, obtained by different methods. 



The full 2MASS designation is 2MASS J19383260-(-4603591. 
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Fig. 1. — Top panel: Our best-fitting light curve model from Bi- 
nary Maker 3.0 (solid black line) plotted on top of a randomly- 
selected cycle in the Kepler light curve (blue points). The sec- 
ondary eclipse observed during this cycle happens to be partic- 
ularly affected by the constructive interference of pulsations and 
chunking (see text). The average photometric accuracy in each 
point is around 0.04%. Bottom panels: Enlarged portions of the 
top panel showing the primary eclipse (left) and secondary eclipse 
(right) profiles. The poor sampling across the eclipses combined 
with chunking and the pulsations leads to an asymmetry in the 
observed profiles that varies from cycle to cycle. 



Here we present an analysis of more than six months 
of short-cadence Kepler photometry taken during QO, 
quarter 5 (Q5) and quarter 6 (Q6) of operations. Our 
primary goal was to measure the secondary time delay. 
Using fits to the eclipse profiles we confirm that there 
is a time delay, which we use in combination with data 
from 0fO to compute masses for the sdB and M dwarf. 
We compare our results with the light-curve modeling 
solution of 010. 

2. OBSERVATIONS 

In May 2009 and from March 2010 to September 2010, 
Kepler observed 1436 orbital cycles of 2M 1938 using 
short-cadence observations. In this operating mode, nine 
6.02-s exposures, each with a readout time of 0.52 s, are 
summed into memory to produce an image every 58.85 
s (92% duty cycle). We downloaded the public QO, Q5 
and Q6 light curves from the Kepler data archive0. We 
converted the time stamps, given as Barycentric Kepler 
Julian Date and accurate to ± .05 s, to Barycentric 
Julian Date. IGilliland etHI (poTol) give additional char- 
acteristics of Kepler short-cadence data. 

3. MEASURING THE R0MER DELAY 

We constructed a model light curve to use as a tem- 
plate for fitting each eclipse, using Binary Maker 3.cQ 

^ http://archive.stsci.edu/kepler 
^ http://www.binarymaker.com/ 
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TABLE 2 


Dates of Secondary Eclipse Minimum 
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with the orbital parameters reported by 01CQ. As shown 
in Figure (TJ the model reproduces the observations rea- 
sonably well. We cross-correlated the primary and sec- 
ondary eclipse template profiles against each observed 
eclipse to determine the times of minima. The model 
was re-sampled step-by-step as it was swept across each 
eclipse to match the sampling of the observations. Since 
finite integration times int roduce asymm etric distortions 
to the eclipse profiles (see iKippind 120 lOl) . we attempted 
to emulate the exposure times by integrating each sam- 
pling point in a similar manner as the Kepler data to 
better match the observed profiles. Before the correla- 
tions were made, each eclipse and its surrounding contin- 
uum were cropped from the light curves and normalized 
by polynomial fits to the continuum. The same fit was 
used for every model-observation pair. We note that we 
do not fit the eclipse depths or durations in this work; 
only times of minima are reported. We also investig ated 
the effects of both orbit al (jShakura fc Postn"ovlll987| ) and 
rotational ()Grooti[2012[ ) Doppler boosting, which asym- 
metrically skew the eclipse profiles, but found only negli- 
gible changes in the fit parameters for the expected levels 
of boosting in this system. 

Tables 1 and 2 (full tables available online) present 
the measured eclipse times, which have an average un- 
certainty of a few seconds. Using the linear orbital 
ephemeris given by 010 as a starting point, we con- 
structed two observed minus calculated (0-C) diagrams, 
one for the primary eclipse and another for the secondary. 
These are shown in Figure [2] and reveal a seven-second 
phase shift between QO and the beginning of Q5; the data 

* 010 do not cite a secondary bolometric albedo; we used 
A2=1.2. 
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Fig. 2. — Top panels: Primary and secondary eclipse O-C di- 
agrams (not on the same scale) from the QO, Q5, and Q6 Kepler 
light curves. Both plots show evidence of a period decrease. The 
dotted line represents 01O's linear ephemeris from QO, while the 
solid line shows our quadratic ephemeris. Bottom panel: Time 
delay of the secondary eclipse with respect to one half period after 
the primary eclipse, which does not significantly change from QO 
to Q5/Q6. The mean AtgE is shown by a solid horizontal line. 



cannot be fitted well with a linear ephemeris. From a 
parabolic fit to all three quarters, we report an ephemeris 
for primary eclipses defined by 

To = 2455 369.422 466 9 ± 0.000 000 5 BJD 
P = 0.125 765 251 ± 0.000 000 002 days 

P = (-1.23 ±0.07) X 10"^" 

The zero-point was chosen to be near the transition time 
between Q5 and Q6. Fits to the secondary eclipse O-C 
diagram (middle panel of Figure[2|) give consistent results 
for the period and its first derivative, which corresponds 
to a nearly 4 ms decrease in the orbital period each year. 
The derived period during Q5/Q6 does not agree with 
the QO result of 010, but this is due to the non-zero 
P, which was unknown at the time. Using our updated 
ephemeris to derive the period during QO, we find the 
same value they do. Orbital deca ys have been obs erved 
for other sdB-l-dM systems fe.g.. lOian et al.l 120121 ) and 
might be explained by tidal dissipation, magnetic brak- 
ing, and gravitational wave emission; the latter effect, 
however, predicts a period change several orders of mag- 
nitude slower than what we observe for 2M 1938. We 
note that a long-period sinusoidal oscillation in the O-C 
diagram could mimic a P; only additional measurements 
over an extended baseline will help discriminate between 
these possibilities. 

As mentioned above, we do not expect the secondary 
eclipse to occur exactly halfway between primary eclipses 
since the mass ratio given by 010 is far from unity. In the 
bottom panel of Figure [2l we plot deviations of the sec- 
ondary eclipse times from 1/2 period after the primary 
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Fig. 3. — Histogram of the secondary eclipse time delays, mea- 
sured with respect to the mid-point between primary eclipses. Data 
represent 1493 pairs of primary and secondary eclipses obtained in 
QO, Q5, and Q6. The Gaussian (solid line) that best fits the dis- 
tribution is centered at 2.06 ± 0.12 s. 



eclipses (hereafter, A^se). Even with the large spread in 
AtsE values, the mean is visibly offset from zero. The off- 
set becomes clearer in a histogram of the measured time 
delays, as shown in Figure [3l A Gaussian centered at 
2.06 ± 0.12 s with a fuU width haff maximum (FWHM) 
of 11 s fits the distribution well. Its centroid agrees with 
the mean of all AtsE values, 1.88 ± 0.13 s. Indeed, the 
secondary eclipse lags behind the mid point between pri- 
mary eclipses, by approximately two seconds. 

Before proceeding, we computed the Fourier trans- 
forms (FTs) of the O-C and AtsE curves to identify and 
analyze any periodic signals that might be present. Only 
Q5 and Q6 were considered in this exercise since the ad- 
dition of QO and the gap it introduces complicates the 
window function. The upper panels of Figure |4] present 
the O-C diagrams next to their amplitude spectra, which 
are plotted out to the Nyquist frequency (3.976 day^^ = 
1/2P). The bottom panels shows the same diagrams for 
the secondary eclipse delay. 

Several strong periodic signals are present in the eclipse 
timings, all of which arise from the rapid pulsations of 
the sdB star and inadequate emulation of the integra- 
tion time of the observations. The latter issue, which 
we refer to as a 'chunking' problem, stems from the fi- 
nite sampling time of each Kepler point, which distorts 
the light curve (see lKipping||2010l ). Even though we at- 
tempted to avoid the effects of chunking by taking inte- 
gration into account in our model template, small mis- 
matches in the way our model and data were binned re- 
main. The sam e problem in accounting f or chunking was 
encountered bv lKipping fc BakosI ()201lD . The effects of 
the sdB pulsations on the eclipse timings are even more 
pronounced. Although none of the 55 pulsational fre- 
quencies detected by 010 has an amplitude greater than 
0.05%, the Kepler data are so precise these small os- 
cillations also introduce noticeable asymmetries into the 
eclipse profiles. The phasing of the distortions changes 
from cycle to cycle, pulling around the best-fitting cen- 
troid values on timescales related to the beating of the 
pulsations with the orbital period. In some cases, such 
as the /=0.26 d^^ signal, the peak-to-peak amplitude 
about the mean approaches several seconds. 

To model the influence of pulsations and chunking on 
the eclipse timings, we constructed a fake, noiseless light 
curve using our binary maker 3.0 template with the 
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Fig. 4. — Top panels: O-C diagrams (left) constructed from primary and secondary eclipse timings in Q5/Q6 and their respective 
Fourier transforms (right). Several periodic signals are visible in the eclipse timings with periods ranging from 0.25-3.7 days Bottom 
panels: Time delay of the secondary eclipse with respect to one half period after the primary eclipse (left) and its FT (right). Some of the 
same periodicities detected in the O-C diagrams are also present here. 
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Fig. 5. — Same as Figure |4l but for a synthesized binary light curve with all 55 pulsation frequencies from 010 added. The pulsations 
and finite sampling skew the eclipse profiles in a periodic fashion that gives rise to apparent oscillations in the arrival times of the eclipses. 
Note the alignment between these signals and those in Figure |4] No noise was added to the model light curve. 
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same sampling and integration as the Q5/Q6 Kepler 
data. We included all of the known pulsations to the 
light curve, using the frequencies and amplitudes given 
in Table Dl of 010. Their phases were defined randomly 
since they were not provided by 010. We forced the sec- 
ondary eclipse to occur exactly halfway between primary 
eclipses to determine whether the pulsations affect our 
measurement of the mean time delay. A second model 
light curve, identical to the first sans pulsations, was also 
constructed. We repeated the entire analysis procedure 
for the synthesized data; Figure [5] shows the resulting 
eclipse measurements and their amplitude spectra. 

The infiuence of pulsations and finite sampling on the 
eclipse timings is drastic; the 0-C diagram FTs reveal a 
large number of periodicities, none of which represents 
a 'real' signal in the data. A comparison of Figures S] 
and [S] shows that each frequency in the observed timings 
corresponds to one of the spurious signals generated by 
adding the finite sampling and the sdB pulsations to our 
model light curve. Most of the signals arise from the stel- 
lar pulsations. Their presence, however, does not offset 
the mean delay measured for the secondary eclipse. The 
best-fitting Gaussian to the distribution of AfgE values, 
shown in the top panel of Figure [Bl is centered at 0.03 ± 
0.11 s, consistent with the modeled offset of zero. Luck- 
ily, the spurious timing oscillations are so rapid compared 
to the run length that we sampled an adequate number 
of cycles during Q5/Q6 that their mean is very nearly 
zero. For this reason, the uncertainties of the pulsation 
phases do not affect our results. If we repeat the above 
analysis using the model without pulsations, the distri- 
bution is significantly narrower, as seen in the bottom 
panel of Figure [6l The best- fitting Gaussian centroid in 
this case is -0.023 ± 0.003 s. Although the pulsations 
don't significantly affect the mean time delay measured, 
they do infiate the distribution width by a factor of 9.7. 

Our light-curve synthesis exercise tells us the centroid 
of the best-fitting Gaussian to the observed AtsE distri- 
bution in Figure [3] refiects the true offset in the system, 
even though the distribution of measurements is infiated 
significantly by the pulsations. Thus, we report a time 
delay of A^se = 2.06±0.12s for the secondary eclipse, as 
measured from the mid point between primary eclipses. 
Note that this offset remains constant from QO to Q5/Q6 
(bottom panel of Figure [2|) , in spite of the seven-second 
phase shift from the period change. 

4. DERIVATION OF THE MASSES 

We can derive the masses of the hot subdwarf (Mgd) 
and cool companion (Mc) fr om the e clipse timings using 
the equations presented by iKaplanI ()2010[) . In a binary 
system with exactly circular orbits and unequal masses, 
secondary eclipses will not occur exactly 1/2 period after 
the primary eclipses due to the extra light-travel time 
(LTT) required. This R0mer delay (AtLTT) is given by 



TTC \q 



(1) 




Time Delay [s] 



Fig. 6. — Same as Figure [S] but for a noiseless, synthesized light 
curve with (top) and without (bottom) pulsations. The model 
was constructed with secondary eclipses occurring exactly halfway 
between primary eclipses. The presence of the sdB's pulsations 
inflates the distribution of eclipse measurements by a factor of 9.7 
but does not change the mean significantly. 

eclipses through an additiv^^ term Ate defined by 

2Pe 



■ cos a; 



(2) 



where P is the orbital period, Ksd the sdB orbital veloc- 
ity, q the mass ratio (Mc/Mgd). Small eccentricities will 
affect the relative timing of the primary and secondary 



where e is the eccentricity and u; the ar gument of pe- 
riapsis (ISternd 119401: IWha[2010l: IKaplanl '201[R). Thus, 
for small eccentricities, the total shift of the secondary 
eclipse with respect to 1/2 orbital period after the pri- 
mary eclipse is 

AisE ^ AtLTT + Ate (3) 

In SjSl we found an average secondary eclipse time de- 
lay of about two seconds with respect to the mid-point 
between primary eclipses. This measurement represents 
AtsE in the above equation, and in order to use it to cal- 
culate the masses, the contribution from the eccentricity 
(Ate) must be identified. Surprisingly, an eccentricity 
of only e = 0.0003 could shift the secondary eclipse by 
the full amount we observe; to separate the R0mer delay 
out from the total measured delayed requires knowing 
the system's eccentricity to a level of precision better 
than this. The radial velocity curve published by 010 
only limits e to 0.02 or smalleiQ. Theoretically, one 

^ Eqn. [T] should also be multiplied by an eccentricity term, but 
this addition changes AtLTT by less than 1% for e < 0.1. We 
ignore it here. 

^ The IKaplanI (120101 ) expression (his Eqn. 6) is missing a factor 
of 2. 

as estimated from an F-test showing that e > 0.02 gives a 
significantly worse fit than a circular-orbit solution. 
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could measure the eccentricity by looking for the apsi- 
dal motion predicted by general relativity and classical 
mechanics, which is straight-forward to calculate. The 
periastron advance (changing uj) induced by these ef- 
fects would give rise to oscillations in the eclipse tim- 
ings, and in particular AtsE, with known periods. The 
presence or absence of such a signal in the data allows 
one to place upper limits on the eccentricity. For typical 
eclipsing sdB-l-dM systems, unfortunately, the expected 
precessional periods are several decades, too long to be 
measured with the current dataset. 

Without knowing the eccentricity to the required pre- 
cision, we continue under the assumption of a circular 
orbit. Combining the observed 2.06-s delay with P and 
Ksd (from 010), we derive a mass ratio (via Eqn. [T]) of 
q = 0.2691 ± 0.0018. This result is independent of the 
orbital inclination angle. Upon assuming a particular 
inclination, we can also solve for the individual masses 
by combining Eq n. [1] with Kepler's Third Law, as done 
bv lKaplanI ()2010l ) in his Eqn. 7. He assumes a perfectly 
edge-on system in deriving this expression, and so it must 
be modified by a multiplicative term (sin for our 
use. If we adopt the light curve modeling result i=69.45 
± 0.02 deg from 010, wc derive masses of Mi =0.372 ± 
0.024 Mq and M2=0.1002 ± 0.0065 M©. Additional or- 
bital parameters calculated from the timing method are 
summarized in Table |3l 

Our derived mass ratio does not agree with the light- 
curve modeling results of 010, from which we infer 
q=0.244 ± 0.008; the disagreement in the masses them- 
selves is even more pronounced. Their mass ratio pre- 
dicts a time delay of 2.35 ± 0.10 s, which is 0.29 ± 0.16 
s longer than our result. This difference (significant at 
roughly 2a) , can easily be explained if the eccentricity is 
as smaU as e = 0.00004 ± 0.00002! We cannot measure 
such a minuscule departure from the circular-orbit case 
using the currently-available data. Although non-zero 
e likely explains the disagreement, shortcomings in the 
light curve modeling might also be at play. Presumably, a 
revised light-curve analysis with q fixed at 0.2691±0.0018 
would result in a slightly different inclination. Unfortu- 
nately, the brief description of the light-curve modeling 
given in 010 does not allow us to predict a revised incli- 
nation with confidence. Our inferred masses are therefore 
provisional. We note, however, that in order to get an 
sdB mass equal to 0.48 M© (01O's derived value) using 
the observed P, A^ltt, and Kgd, the inclination would 
have to be 66.7 deg, in disagreement with the light curve 
modeling result. This discrepancy might be explained 
by a number of uncertainties associated with the light 
curve solution, including how to accurately model the M 
dwarf albedo and the Kepler bandpass. Their solution 
also depends greatly on the quoted value of log g. Sig- 
nificant changes larger than their 1-cr error (0.009 dex) 
might come about from non-LTE modeling (versus their 
LTE solution), unaccounted-for rotational broadening of 
the H Balmer lines, and incomplete removal of the orbital 
velocities upon summation of the spectra for fitting. In 
light of some of these effects, 010 themselves stress that 
their model parameters are subject to systematic errors 
that might be larger than the statistical errors they cite. 



5. CONCLUSIONS 

Using three months of short-cadence Kepler photom- 
etry, we have shown that primary and secondary eclipse 
timings can help constrain the component masses in 
sdB-|-dM ecl ipsing binaries using the technique described 
by iKaplanI ()2010i) . The secondary eclipses of 2M 1938 
arrive nearly two seconds after the halfway point be- 
tween primary eclipses. Assuming a circular orbit, we are 
able to derive the mass ratio and individual component 
masses from this delay, the orbital period, the sdB veloc- 
ity, and the inclination angle. Our total system mass and 
the mass ratio disagree with the results of 010. However, 
this difference may be reconciled if the system has an ec- 
centricity no less than '--^0. 00004, which is too small to be 
measured using the currently-available d ata. This work, 
and th e concurrent study of KOI- 74 by iBloemen et al.l 
(|2012D , represent the first detections of secondary eclipse 
R0mer delays in compact binaries. The results of both 
studies, however, are fundamentally inconclusive at this 
time; until the eccentricities in these systems are deter- 
mined to a precision of ~ 10~^, it is impossible to sep- 
arate out the R0mer effect contribution from the total 
observed time delay. An eccentricity-induced time delay 
can easily and naturally explain the apparent disagree- 
ment between the R0mer delay and binary light curve 
modehng solutions for both 2M 1938-^4603 and KOI- 74. 

If the binary light curve modeling solution and derived 
mass ratio from 010 are accurate, the time delay we 
measure implies 2M 1938 has a non-circular orbit with 
eccentricity no smaller than e ~ 0.00004. Close sdB-l-dM 
binaries are generally expected to have circular orbits, 
as they are likely the products of CE evolution. Even if 
some level of eccentricity remains after their formation, 
their circularization timescales were thought to be short 
enough for the orbits to circ ularize within t he lifetime o f 
the sdB, around 100 Myr (jTassoull [19881 [Zah^ [T971 . 
Several phenomena might explain a small eccentricity 
in 2M 1938' s orbit. Soine rec ent studies of sdB+dM 
systems (e.g.. fPablo et al1l2012t ) show the synchroniza- 
tion timescales might be longer than the EHB lifetime. 
The system also might have been born (following the 
CE phase) with the M dwarf's spin misaligned with the 
orbit, since in the wider pre-CE binary there was no rea- 
son for it to become aligned. Spin-orbit misalignment 
has been reported for the young eclipsing binary DI Her 
(jAlbrecht et al.l [2"009f ). The resultant torques will give 
rise to a time-varying eccentricity, apsidal motion, pre- 
cession of the orbital plane, and spin precession. (In this 
case the misaligned figure of the M dwarf would probably 
vitiate some of the assumptions of the light-curve model 
used by 010.) A third body orbiting outside the 3.0 h 
binary would perturb any initially circular orbit, result- 
ing in orbital eccentricity. If the third-body orbit were 
not coplanar with the inner binary, one would expect to 
observe precession of the inner binary's orbital plane as 
well. Any eccentricity will be damped out, at the ex- 
pense of orbital energy, so another prediction is that the 
orbital period should be decreasing with time (as ob- 
served, although a longer timebase is needed to be sure 
that this is a truly secul ar phenomenon, rather than pe- 
riodic or quasi-periodic) . lEggletoiJ ()2006D discusses these 
processes and their associated amplitudes and timescales 
extensively. Random convective motions in the M dwarf 
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TABLE 3 
System Parameters 



Par am 


Value 


Error 


Unit 


Comments 


P 


3.018 366 023 


± 0.000 000 048 


hrs 


measured from Kepler Q0+Q5+Q6 primary minima at start of Q6 


P 


-1.23 


± 0.07 


10-10 s s-i 


measured from Kepler Q0+Q5+Q6 primary minima 




65.7 


± 0.6 


km s""'" 


measured bv 0stensen et al. ('2010^ 


i 


69.45 


± 0.02 


dcg 


caleulated bv 0stensen et al. f2010) 


AtsE 


2.06 


± 0.12 


s 


measured from Kepler Q0+Q5+Q6 (from Gaussian fit to distribution) 


9 


0.2691" 


± 0.0018 




calculated from A^se, P, Kg^ 




0.372° 


± 0.024 


Mo 


calculated from AfsE, P, Kgdi sin i 


Me 


0.1001° 


± 0.0065 


Mq 


calculated from AisEi P^ K^di sin i 


a 


0.823° 


± 0.015 


Rq 


calculated from P, K^d, M^d+Mc 




0.196° 


± 0.049 


Ro 


estimated using log g and M^d 




244.2° 


± 2.8 


km s "'" 


estimated using q and iCgd 



"assumes a circular orbit. 



could also lead to an eccentric orbit, since they induce 
flu ctuatio ns in the exterior g ravitational field (Phinney 
[T992t iLanza fc Rodon6l[200l . 

Due to strong irradiation, we expect some level of H- 
alpha emission from the M dwarf. If such a feature 
is detected and the eccentricity of the system is even- 
tually determined to adequate precision, 2M 1938 will 
be unique among binaries in that the individual masses 
can be computed using at least four different techniques: 
double-lined spectroscopic binary analyses, asteroseis- 
mology, eclipsing binary light curve modeling, and eclipse 
timing monitoring. Comparisons of these results will 
help constrain the underlying physics in the light-curve 
modeling and asteroseismic codes. 
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